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Abstract 


A theoretical investigation is undertaken into the creep 
instability of two common shells of revolution, the infinitely 
long cylinder and the complete sphere. The linear and non- 
linear visco-~elastic material idealizations are used in the 
various theoretical procedures. The critical buckling time 
is greatly influenced by the parameter a/h ratio, magnitude 
of initial imperfections and material selection. The spherical 
Shell has a lifetime longer than the cylindrical shell having 
the same a/h ratio, initial imperfections and material con- 
stants, while the critical mode numbers are the same. 
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NOTATION 


W transverse displacement, positive outwards. 
u displacement inn direction. 
1 

V displacement in ue al vectron. 
M ,M ae moments per unit length as defined in Figure 2. 

2 
‘Oe forces per unit length as defined in Figure 2. 
e ,€ ,€_ ,e total strains. 

meee 2 21 i12 
a external loads per unit area as defined in Figure 2. 
ere = orthogonal coordinates as defined in Figure l. 
O70 transverse shear forces per unit length as de- 

} 2 fined in Figure 2. 
Rae Principle radii of curvature of a shell and 

1 2 radii of curvature of the coordinate curves 

Nn 7 - 
LZ 
a cylindrical shell mean radius. 
an out-of-roundness parameter defined by equa- 
elon, (4. 1o)e 

h cylindrical shell mean thickness. 
S arc length. 
A,B defined in equations (2.5 b-c). 
B distance measured along a normal from the 


middle surface of a shell. 


E Younos Moadulussor Ellasticity- 

C,K material viscous coefficients. 

14 creep exponent defined by equation (3.4). 

R spherical shell mean radius 

(2 out-of-roundness parameter of spherical shell. 





NOTATION (cont.) 


spherical shell mean thickness. 

amplification function (cylinder). 

time. 

CYrLtELcalmtime - 

amplification function (sphere). 

maximum initial imperfection; cylinder/sphere. 
harmonic in circumferential (8) direction. 
Deomiomienan cireumpolar (Od) direction, 

lines of principle curvature. 

change of curvature. 


half total change of curvature in twist. 





ae I 
principle curvatures = a : 


| eer 


f 


membrane strains. 

shear strain. 

DOlSSOn S Katio- 

effective stress defined in equation (3.5b). 
4 4 

3 Gn” OT 3 GH. 


circumferential and axial coordinates in 
eylindyvical sheild. 


circumpolar and circumferential coordinates in 
spherical shell. 





= kh or = ki. 
Ch3 cH 
joo oa 
Kh3 KH3 
er 





NOTATION (cont.) 


-aaw/E (time constant linear viscous cylinder). 


-aQw/T (time constant non-linear viscous 
eviinder)-= 


-3Row/2E (time constant linear viscous sphere). 


= oR / 2G (time constant non-linear viscous 
sphere). 


xt (non-dimensional time linear viscous cylinder). 


Xt (non-dimensional time non-linear viscous 
cylinder). 


dt (non-dimensional time linear viscous sphere). 


At (non-dimenSional time non-linear viscous 
sphere). 


n(n + 1) 
defined by equation (5.12c). 
O*/0< a. 
0 
Stress voevlaLor. 


initial value of ( ). 


oa) 
di ° 





J) 
dt : 





dominant part of ( ). 


perturbablon spare of ()))- 
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1. INTRODUCTION 


ile Preface 

The problems of creep deformation and buckling of structural 
members while often addressed in aerospace and chemical engineering 
have, until recently, been aphonic design considerations in ocean 
engineering. The rationale being that the ocean, itself, offers an 
environment very unaccommodating to creep with its relatively low tem- 
peratures. Further, the general use of creep resistive materials for 
relatively short periods of loading made creep insignificant in the de- 
Sign of ocean shell structures. Recently, the escalating démand man 
and his society have placed on the ocean in pursuit of its understand- 
ing and resources has resulted in a rapid advancement in ocean technology. 
This new technology has generated an increased interest in the creep of 
shell structures. This manifestation of interest has been stimulated 
by many events, the greater use of creep prone materials, aluminum and 
titanium, the longer duration structures are required to withstand com- 
pelling hydrostatic pressures, and the need to predict volumetric changes 
due to creep of shells which have stringent volumetric design constraints, 
such as the neutrally buoyant and compressible, free floating, current 
measuring SOFAR floats of the Woods Hole Oceanographic Institution. 

In the past, the ocean engineer was perplexed with a two part 
dilemma when asked to predict the creep performance of his design, first, 
find ‘an acceptable analytical model for creep deformation of shells and 


second, locate accurate material constants for low temperature creep. 


o= 


The purpose of this paper is to solve the first dilemma of 
developing a simple analytic technique to predict creep deformation of 
two common shells of revolution, the cylinder and complete sphere, suit- 
able to oceanic applications and which will be useful in preliminary 
stages of design and may be adequate for final design. Finally, to dis- 
cover if creep buckling is a phenomenon that needs to be actively con- 
Sidered in the design of long term ocean shell structures. 

1.2 Historical Review 

A literary search of the existing literature reveals that few 
published works have been transcribed which are directly applicable to 
conditions experienced in the ocean. Most papers deal with colums, 
plates or cylindrical shells loaded in axial compression. A number of 
investigators have attempted solutions for the creep behavior of thin 
cylindrical shells loaded ager uniform extermal pressure. To the best 
knowledge of this author, no paper has been published dealing with the 
creep performance of a complete spherical shell loaded under external 
pressure. 

The author's dealing with cylindrical shells use a variety 
of simplifications, but do have an assemblage of assumptions common to 
all their analyses. These similar assumptions include: 

(a) thin shell assumption, so that the radial stress compo- 
nent is negligible in comparison to the in-plane stress components 

(b) only the effect of secondary creep, according to Odquist's 
law (1937) are considered 

(c) buckling is not defined in the bifurcation sense, but 


rather as a state of progressive accelerated growth rate of the displacement 
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field, creating unfavorable strain distribution and leading to full 
Serpuccural failure. 

(a) initial imperfections of a shell cause eccentricity of 
loading, thus serving as a lever arm for compressive in-plane forces, 
which in turn increases the bending moment, which in turn, increases the 
rate of curvature change 

(e) an assumed solution of an infinite trigonometric series, 
either a complete or partial Fourier series. 

Hoff, et al (1959) used a sandwich model, in which the outer 
skins support membrane and bending stresses and the core supports shear 
without deformation. The sandwich model avoids soup ion tiene in the 
superposition of compressive and bending stresses due to the non- 
linearity of the creep law. Further, a quasi-elliptical initial imper- 
fection was posure Samuelson (1970) developed a numerical solution 
for non-uniform external loads. Sankaranarayan (1969) generated a solu- 
tion for creep deformation of cylindrical shells under combined lateral 
and axial pressures for a material that adheres to the Tresca criterion. 
Berman, et al (1974) presented a numerical analysis for a cylindrical 
shell infinitely long, so that, a state of generalized plane strain exists. 
As did Hoff, Berman assumed that the initial imperfection is quasi-elliptical. 

Penny and Marriott (1969) have carried out the only known, but 
rather special analysis of creep deformation of a spherical shell. They 
explored snap through of a spherical boss point loaded. 


All analyses emphasize the fact that initial imperfections 


play a most important role in creep buckling. That, buckling stems 
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from the accelerated growth of these initial imperfections in an otherwise 
uniform initial displacement field. 
ieoeewethod of Analysis 

The conviction that creep buckling stems from the growth of 
small imperfections in the otherwise uniform initial displacement field 
suggests the use of a perturbation method of solution as first developed 
by Abrahamson and Goodier (1962) for dynamic plastic problems. The 
method assumes an axisymmetric dominant solution in which buckling does 
not occur and a perturbation solution in which the influence of bending 
is significantly greater than the in-plane forces thereby giving rise 
to inextensible buckling. 

Due to the axisymmetry of the dominant displacement field of 
the shell of revolution, only the transverse displacement field consists 
of a dominant and perturbation part, while the in-plane displacements 


have only a perturbation part so that 


wientioe Ee) =awilt) ow! (ny mort) (1.1a) 
i) eae 1 2 
eG peerage ent ilies 9,2) (1. 1b) 
ee 1 2 
and 
Geet 7 C) = verre, Mav ©) Cie) 
i 1 22 


respectively, with perturbation parts delineated by primes and dominant 
parts by bars. 
In turn, to remain consistent, the shear strain and twisting 


curvature change have only perturbation parts while the membrane strains 


and bending curvature changes have both dominant and perturbation parts. 


ee 


Similarly, the shear forces and twisting moments have just perturba- 
tion parts while normal forces and bending moments have both dominant 


and perturbation parts 


N@m,n, ef =N (t) +N' mM, 7, t) (ee) 
‘eel l ie 1 2 
Meh: Tee = N Ceyeen (ne, 17 ©) (1.2b) 
De Ae 2 Dn 2 
Cech) pei) = NO Gemeente) = NiCr ee) (1226) 
nay. ) Wapk oF Di a ae 
Mi(n,n, t) =M (t) +M' (1,7, t) (Guan 
ie | 2 1 ioe 2 
Min,n, t) = M(t) +M' (n,n, t) Gees 
Z. l Zz WA Pe i] 2 
and 
ea te c= My) > No, ee Mines Mee 2b) (2) 
i 1 2 io 2 a Gee 


The material comprising the shell is isotropic, incompressible 
and homogeneous. The initial elastic response is considered to be infini- 
tesimal. No initial shape of imperfections is inferred. Further, the 
assemblage of assumptions common to past analysis, (a) - (e) on page ll 
is assumed. 

The analysis will present non-linear equilibrium, strain- 
curvature relations and constitutive equations for linear and non- 
linear visco-elastic materials for an arbitrary shaped shell. The 
reduction of these relations will be carried out for the infinite cylinder 


and the complete sphere. 
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2. BASIC RELATIONS FOR THIN SHELLS 


fie introduction 
The succeeding relations were developed by Walters and Jones 
(1972). They were derived with the aid of differential geometry as 
discussed in Kraus (1967) and are valid for any thin arbitrary shell 
which has lines of principle curvature aligned with the parametric 
lines. 
2.2 Shell Coordinates - Differential Geometry 
The coordinates of point P on the middle surface of a shell, 
the reference surface, are defined in Figure l. r is the position 
vector which is a function of parametric lines ne es AS Es the unit 
normal vector. The poSition vector can be expressed as 
AN 


ra xd + xQe xk 
1 Z 3 


(2.1) 


a7 ae ' ; 
tee 7, and kK are the unit vectors in the x , x , x directions, 
i Z 3 
respectively. 
As demonstrated in Kraus (1967), the magnitude of the prin- 


ciple curvatures are the roots of the following quadratic equation 


H2“2 — (EN-2MF-GL)K + (LN-M2) = 0 (2) 
where 
E=f, ‘Yr, (25 3a) 
1 1 
F=f, f, (2.35) 
1 2 


(equal to zero when n is orthogonal to n ) 
} Z 


=15= 





G.= 1,8, (2. 3c) 


je D) 
1. 
H = (EG-F2) ” (2.3) 
Pe 
ic. = reise (22.2€) 
el 
fee 
Mo ne (2 ous) 
1 


(equal to zero when n= and as coincide with the lines of curvature) 
Il 


“N 
N=ner, (2.34) 
with _- — 


= 7 (224) 


Equations 5{a-d) are referred to as the first fundamental magnitudes 
while equations 5(e-g) are called the second fundamental magnitudes. 
ile LF and is are orthogonal, then F equals zero, the 


elemental area can be expressed as 


ds* = a*a2n + B2da2n (2.5a) 
1 1 
where 
= 45 
A=) (2.555) 
and 
35 
B = (G) (2 soe) 


2.3 Strain-Curvature Relations 
‘ A tensor formulation and the usual assumptions of shell 
theory, (eg. Kraus (1967)) have been employed by Walters and Jones (1972) 


to show that the strains in the thin shell may be written 


=16= 





e = £ eae (2.6a) 
ot l | 


e =e + 4K (2.6) 
a2 2 2 


S ae 


(2.6c) 
12 2 


lI 
= 
=) 
+ 
NO 
WN 
w 


where the membrane strains are 








u, Ay Vv (w, )? 
l 2 “6 
is = el + AB a “a (2.7a) 
l 1 22 
Yi ee (Wiese 
Z l “7 
S = + = ’ (2. 7b) 
y) B AB - op2 
and 
ee 
A,u B,V 
= —(— + —(— ae 25 
La Ba), a'B , a AB (2276) 
Z 





W, A, Wy 
l.u l 2 Vv 2 
oe eee (a 2.74 
1 aR A - AB o B ( ) 
I) l 2 
W, ee 7s 
er, 4 2 u 
SS SUR ~ Bp? * HBR” a? (2. 7e) 
2 
and 
wos Wy 
i Boy 
Pa = —(——— — pede ee, (2275) 
BAR B 
12 Pa pe 2 eee 


Note that if is the total shear strain while co is one-half the total 


twist curvature change 


5 


The Love-Kirchhoff assumption has been invoked in the 


derivation of these equations. 
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2.4 Equilibrium Equations 
The typical shell differential element illustrated in Figure 2 
will have a rate of internal energy dissipation. 


W. ffi sudan y emt 
a ie 22 12 0 11 


+MkK +2M kc. (Idas (2.8) 
22 12 12 


where s is the total undeformed area of the middle surface of the 
arbitrary thin shell. Now, equating the rate of internal energy dis- 
Sipation to the rate of external work (i.e., using the principle of 


virtual velocities), a consistent set of equations of equilibrium is 














obtained. 
ABp - NB, + (N A), + (NB), +N A, 
l 2° See eee 2 bl ae 2 
(MB), MB, (M A), M A, 
] ] Vn | 1:2 Zz 2) 2 
1 ee +---_---- = O (2.9a) 
R R R R 
] ] ] i 
ABp - NA, + (N B), + (N A), +N B, 
Zz | eZ 2 ] 2 Z 2 1 
(M A), MA, (M B), M B, 
2 ee 2 Pe ] se 
+ — - Be 
DP R a O (27, 9b) 
2 2: 2 2 
and 
N AB N AB N Bw, N Aw, 
AB - - + + ; 
ee x R A oF B 2 
2 


(MB), 
yt 


Fein wee, “Ne i ee 
2 oer) tie 2 





12 (2 ioe 
+ + 
B 9 ¢ A “7 oe B 15; 
MA, MB, 
5a 
+ 2M 2 | ae = 6 (2.9c) 
In B /'9 A }*y 


Inertia terms are disregarded since creep is assumed to be a very pro- 


tracted deformation process 


3. CONSTITUTIVE RELATIONS 


3.1 Linear Visco-elastic Material 

The appropriate model to simulate secondary creep during the 
application of a step function load is the Maxwell model illustrated 
in Figure 3a. The strain at all times is the sum of the elastic strain 
component, and viscous strain component. The stress in the viscous 
and elastic elements are always equal. Thus, from the equation of 
equilibrium the following expression for strain rate is found, 
(eg., Drucker 1967) 


4 (3.1) 


(D e 
I 
t]a. 


g 
G 
During application of the step function load, the motion of 
the dashpot is negligible. Hence, the initial strain is the elastic 
stretch of the spring. As time procedes, the model will continue to 
@meeep Steadily at the rate of 0/C as shown in Figure 3b. 
The relation for the linear Maxwell model can be written 


in tensor notation as, (eg., Drucker 1967) 


Qe 
— 


; i * (l= 2U) kk eee 
Pe = 44 
LJ B14 13 3 : 2 eee 
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where 


O = 6) + O + O 
KK es 22 co 
) 
kk 
S =n6) - ——4§ 


alg 13) 32 Lj 


I 


ie kronecker delta 
and E, C are material elastic and viscous constants 


respectively. Since superposition of hydrostatic stresses is assumed 


to cause no further straining,equation (3.2a) reduces to 


je 
-- 
<= 


(32/25) 


Ul 
| 
4- 
Nw 
a 
Qik 


Inverting equation (3.2b) and performing the proper intergration reveals 
the linear viscous constitutive equations, after the initial elastic 


response, that is, disregarding the elastic behavior, for a thin shell as 


N =a(e +=) (3.3a) 
y 1 NS) 
° 1 e 
N =a(e +=>€) (3. 3b) 
2 oe 
a Ld 
N = N = — Cesc) 
11922 Za) 3 Yo 
* 1 s 
M = E(K +=€K ) (3,30) 
l en) 
e iG ° 
M = E(K + —K ) (2 3e) 
g) y) 2 
and 
ie. a 12 
where 
4 
Qa= 3ch 
d ; 
a a ene 
‘ 9 
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3.2 Non-Linear Visco-Elastic Material 
Secondary creep for the non-linear visco-elastic material 


can be, as in the linear case, represented by a Maxwell model where 
(324) 


o.M 
(z) 


De 
Il 

I] Q- 
++ 


This model does not simulate primary creep and has a poor 


response to variable stress conditions. 
Assuming that the creep strains are dependent on loading 


history, the Prandtl-Reuss Incremental relations can be used for the 


non-elastic part 
woe a) 


where 
o* and e* are effective stress and strain based on Von 


fees criterion 
)2 


i 


where 
L p) yh 
o* =— |(o0 -o )*+ + (0 -o )* + (GO -O 
vo 11 22 2 33 33 
+ 6(T Zar 447 2) |" Coo) 
12 23 3] 
ex ee |‘ -e )2 + (e -e 2 + (e -e 2 
/> ieee 27a el 
I 
+ 6(y 24y 44y 2)|* (3.5c) 
LZ 23 ou 


Norton's Law for secondary creep under a uniform constant load 


(3.6) 


1s expressed as 
all 


Gls = z at. 
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Hence, one can substitute (3.6) into equation (3.5) and find 


gS yoL 
0 


° 


ae 
+3 





3 
=> (32) 


K 43 
which is Odquist's Law for secondary creep, (eg., Odquist (1956)). 

By asSuming that the effective stress remains constant across 
the shell thickness, equation (3.7) can be inverted and intergrated to 


yield the non-linear constitutive relations while neglecting the 


elastic effects (Penney and Mariott 1971) 


* = e e 
N = Qo (1 we + LZ: ) (3.8a) 
1 Vig 2 ae 
* law e e 
hee. ic = ee + ae (33ab)) 
g 208 2 al 
* (1-M) 
N =N = Lees (3. 8c) 
12 Za 3 0 
* — e e 
Me = co eC ee) (3.8d) 
1 i €24092 
* — e e 
M = 2n6) (1 Mig + aE ) (3.8e) 
D 02 eal 
ke ae e 
M = =—M = 6 (1 M) (3.8£) 
1 zy) 12 
where 
4 
2? = — kh 
3 
as 


Note, that if M equals one, the linear case, these constitutive 


elactions reduce to the linear relations, equation (3.3). 
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4. INFINITELY LONG CYLINDRICAL SHELL 


fee Introduction 

The objective of the following section is to reduce the 
general equations for an arbitrary shaped shell to the particular case 
Sema infinitely long, thin cylinder. The perturbation solution for 
transverse deflection is assumed to be a Fourier cosine series 


OO 


vo A(t) cos ng (4.1) 


n=0 


@22 Basic Cylindrical Shell Equations 
4.2.1 Differential Geometry 
Defining x as n and 98 as n , as illustrated in Figure 4, the 
l 2 


position vector equals 


— . i as ae 
Yr = Xi + asinOj + acos6k : (4.2a) 


The principle curvature magnitudes are 


* 


ene (4.2b) 
x 


* 


ee lfa (4.2c) 


emo the differential area is 


ds? = dx* + a%a04 (4.2a) 
where, 
A= 1 
and 
B= a 


fmeze2 Strain Curvature Relations 


Equations (2.7 a-f) reduce to 


-~2 3— 








Ee Bey 
Ww, 2 
C= alee + wt : ) 
0 a 8 DBE 
- tte ; . Wr Wr 
an Vy a 
K = -w, 
x xX 
il 
So aed een clog) 
a 
and 
1 
ay ea ies 7 SU ay 
The time derivatives of (4.3a-f) are 
nae Ur, + Wr Wry 
° 1 e ° all e 
Ey ae (Vy, +wt 3 WroWro) 
ee ee ee 
is a <6 x a, ae 0 
< = =w, 
x xX 
« in e ° 
mE I 5 (Vio — Wroo) 
a 
and 
e Hi ° ° 
TAS a ; 2 
x9 a SY ae ) 


4.2.3 Equilibrium Equations 


Equations (2.9 a-c) yield in-plane equations 


+ + alN 
OP One ON x) x 


and 


a Fe 


Ww, 


) 
x 


(4.3a) 


(4. 3b) 


(4.3c) 


(4.3d) 


(4.3e) 


(4.3£) 


(Cie4a 


(4. 4b) 


(Ane) 


(4. 4d) 


(4.4e) 


(4.4f) 


(45a) 
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+ aN + += = 4.5b 
AP 5 oe se No 6 = M + M 0 ( ) 
and transverse equation 


a -~ N. + a(N w + N W 
1 SW os ji) 


6 


+ diy + N w,.) + 2N 
A O° wee ceo 7o x 


0 (4. 56) 
4.3 Buckling Solution 

For the axisymmetric loading of an infinitely long cylinder, 
the following assumptions are stipulated: 

(1) loading is axisymmetric, thus the dominant solution has 


no 8 dependence and ie ha 0 


N 
sue : 6 
(2) plane strain conditions exist, N.. = Sy 
(3) transverse displacement, forces and moments do not vary 
Significantly with the longitudinal, x, direction, hence both the dominant 


and perturbation solutions are independent of x. 


The equilibrium equations simplify to dominant equation of 


» 1 on| 


Se = 0 (4.6a) 


and perturbation equation of 


1 t cs 
Now ry + M 6,06 = O (4.6b) 


Mmendominant rates o£f membrane strains are 


e ean (4.7a) 





(4. 7b) 


- 8 
0 a 
and 
y =0 (4.7c) 
0 
and the rate of curvature change become 
a= 0 K, = 0 a= 0 


a= 
Assuming inextensible buckling the perturbation membrane strains 


are zero and the curvature changes are Simplified to 


K' = 0 (4. 7d) 
x 
1 mt, tn, eee oa aA a 
K' (w W ro9) : (4.7e) 
a 
K' . = 0 (4.7£) 
x0 


4.3.1 Linear Visco-Elastic Case 
Substituting the constitutive relations (3.3 a-f) into 


(4.6a) reveals a dominant transverse displacement rate of 
e a Pp 
= (4.8) 





w= 
The perturbation solution (4.6b), with the assumption of 


inextensible buckling (€' = e'. = «' = 0), becomes 
x S) x8 


“8008 


OWW' ya _ - E/alw',., + w' (4.9) 


) = 0 


(4.1), into equation 


Substituting the assumed solution for w', 


(4.9) yields a first order, ordinary differential equation of constant 


@eerricients. 
Vie fa(n2 21) A =O (4.10a) 
n n 


% 


where 
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AL = amplification function 


xe, (4.10b) 
E 
Thus: 
c 
A = C_ exp eS (4.11) 
n n 2 
(ae) 
where 
< = an arbitrary constant dependent on the 
initial imperfections. 
Defining 
“™ 
T = xt 
- x 
one can plot A /C versus T to discover the critical mode in, 2s 
gle sg @. CEE 


Figure 5 reveals that the elliptical mode, n equal to 2, is the critical 
mode. The growth of the critical amplification function is displayed 
mieragure 6. 
mmemeeenon—-Linear Visco-Elastic Case 

Combining the equilibrium equation (4.6 a-b) and the non- 
linear constitutive equations (3.8 a-f) generates a dominant transverse 


isco lacement rate of 
P y) ‘ (M—1) 
cs a Po 


2 7 ne 
Ww 5 (4212) 


The perturbation equation is transformed to 


Oirwy t ee "4 rae = E 3 
OWT gg” GAN W nae Y eaans) ae ie 
Note, if creep exponent M equals one, equation (4.13) reduces to the 
linear equation (4.8). 

Substituting the assumed solution for w', equation (4.1), into 


equation (4.13) yields a first order ordinary differential equation of 


gg 





constant coefficients. 





ee = 20 (4.14a) 
n n 
where 
A, = amplification function 
and ; 
eee (4.14b) 
(a 
hence 
+ 
AT =: C. exp aa (4.15) 
(n?-1) 


where Co is an arbitrary constant dependent on initial imperfections 


and 


As in the linear case, the elliptical mode is critical. 

In the formulation of the analysis, it was assumed that the 
effective stress remained constant, however, as the perturbation solu- 
tion of transverse displacement grows, the effective stress increases. 

To overcome this problem one can divide the deformation process into a 
finite number of time segments in which the effective stress is considered 
constant. This operation generates a piecewise approximation of the 
transverse deflection. 

Plane strain conditions allows the effective stress to be ex- 


fjressed as 


o* = 0.866 (5, Age) (4.16) 


6 


where 
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Oo. = (Aggy) 


which reduces to 





oa ease), 
6 2 
a 
Defining 
i o* 
ae BG 
0 
the perturbation solution becomes 
Te (we-1) 
wi = Ce exp (> Ww eeCOseu (4.18a) 
and 
; (ee a 
es ee eesans SMa 5 
WwW 5 = Ww i (4.18b) 
The effective stress becomes 
P,4 (M-1) a "oa (Mel) 
o* = .866 aa 1+6W (") a*exp =e (4.19) 
where 


a* c/a (initial out of roundness pa- 
rameter << 1). 


The anitial GfEecEive stress is 


T =X Te 
C 0 
a = initial value of X 
Da 
ra a 
o* = ,866 —— (1 +6 (=) a*) (4320) 
0 h ia] 


-29- 





The duration of the time segment is dictated by the growth of 
effective stress. When the effective stress increases by a given per- 
cent, i.e. ten percent (WW = 1.1), a new time segment using the new value 
meeo* 1S initiated. 

As an example, taking (a/h) equal to 10 and a* equal to .005, 
the amplification function's accelerated growth causes buckling at ee 
Semel tO 2.59, as graphically portrayed in Figure 7. Inspection of 
Figure 7 shows that the amplification function growth is very slow. It 
tenwonly over the final 10 percent of the total lifetime that the ampli- 


fication function becomes substantial, indicating catastrophic buckling. 


» 
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>. _ COMP LEDER esr NE RECA Sul ine 


fe Introduction 
The succeeding chapter presents the creep response for 
the complete spherical shell. The perturbation solution for trans- 


merce deflection is assumed to be 


le @] 


n 
w' = , (eq (01P, (ose + Baa (#) Fe (cos) cos) 


= m=1 (5.1) 


where P and Pe are Legendre and associated Legendre polynomials of 
degree n and order m. 
pee Basic Spherical Shell Equations 
5.2.1 Differential Geometry 
Stipulating $ as oF and 6 as noe as exhibited in Figure 8, 


defines the position vector equal to 


r= R(singo SAGES: cos sine 4 - sind“k) NSieza) 


the principle curvature magnatudes equal to 


K Sy Palo. 
and 
i: il 


and the differential area equal to 


ds? = R2do2 + R@sin?oa62 (5.2a) 


where, 


—-3]- 





and 


5) Cee 


B 


R sing 


Serain Curvature Relations 


Equations (2.7 a-f) simolify to membrane strains 


and 


and curvature changes 


and 


2K 


90 





dhe time derivatives of (5.3 a-f) are 
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Zz 
e se ° csc~o ° 
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(w, 


6” "0 


(5234) 


(S25) 


(5.30) 


(S36) 


(5.3e) 


moe) 


(24a) 


(S5.4b) 


+ ving (Be 
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cS = ren ee 
$6 2 'd 'O6 
e ° ° t ® 
Ko = esc (Vi, - Wr gcc) + cot (u-w, 
R2 R2 
ana 
oe Soe) eee eee 
Te nae ad ee 
- Ww, escd + Qw, cscocotd) 


66 6g 


5.2.3 Equilibrium Equations 


Equations (2.9 a-c) simplify into in-plane equations 


R —- N@cotodot+ N csco + N + N,cot 
pe at coor Ne waeee’ 
+ 2 +M td - M,cotd + M csco) 
R ‘ 6,6 mig Q 66,6 
+ + + N 
Rpé@ Neg, ¢ Seon csco 6,6 
1 
+-+—— 
= (cScOM, 4 + ee + BAY 


and transverse equation 


Zz 
R Po +N Swi Cola +w,,,) +N Wr gp cscd 


> > ) 


- oy) 
eat Ny) 2cscoNn 


+ WN WwW, 


o6°'66  $,6°"6 


+ csc%d + csco(N ) 


Dyer cular BO NerCn | EHORG! 0 


+ M T._ZeCocon 


a Z 
rey Ms + csc~oM 


o,¢ 0,09 


+ 2M cotodcscd + 2M 


66,0 esco - M 


cotd¢ 


99, $9 0,9 
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(5.4d) 


(5.4e) 


(5.4£) 


(Ss. 5a) 


(S,5b) 


(S25e) 





wee Buckling Solution 
In the analysis of the complete sphere, it is assumed the 
dominant solution has point symmetry, loaded in a uniform pressure 


field. As a consequence, there is no dominant 9 or ¢@ dependence, 


= Q, M,=M and N. =N.. The dominant transverse 


oh ? G" 0 > 


equilibrium equation Simplifies to 


p@=pd=0, M, 


Rp, - 2, = 0 (5.6a) 


The in-plane perturbation equilibrium equations reduce to 


M' TeGoroMe 


ane - cotoM’ |, + cscoM' = 0 (5.65) 


> 90,0 


and 


CSCOM 6, 6 + Mae + 2cotdM $6 = QO (57 6c) 


The transverse perturbation equilibrium equation becomes 


NeW reg + Com. + csc*$w', 94) + as 
’ = ' 2 ' ' 
+ 2cotdM Me M 4 + csc ~oM 6,00 + 2cotdcscoM 66,0 
+ Se eerie - SSE + M" = 0 (5.6a) 


The dominant rates of membrane strains and curvatures 
changes are 
a ae = ; Ee = 0 (5.7a-c) 


and 
a = ee , K 40 = 0 (5.7da-£) 


Assuming inextensible buckling the perturbation rates of membrane strains 


are zero and the rates of curvature changes can be expressed as 
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i (w! ure, (5. 8a) 
— a a + csc* bw" a, + cotgw', 4) (5. 8b) 

and 
2K" a = =~ (cscgcotdw",p - cscow'y 46) (5. 8c) 


2 
R 
Tf one examines the perturbation expressions derived for cur- 


vature changes when inextensible buckling is assumed (equations 5.8a - 
5.8c), one discovers that they are identical to the curvature change 
expressions derived by Jones and Ahn (1974). By comparing the formu- 
lation of the transverse equilibrium equations of Walters and Jones and 


a) 


Ahn, it is found that Walters and Jones assumed the term to be neg~ 
R2 | 

ligible in the curvature change expressions, while Jones and Ahn include 
it. The curvature change expressions are assumed to be a function of 
u, v, and w by Walters and Jones, while Jones and Ahn assume curvature 
change expressions are just a function of w. If, however, inextensible 
buckling is assumed, one can express the in-plane displacements, u and 
v as functions of the transverse displacement w only, and arrive at 
the expressions derived by Jones and Ahn. Thus, it is reasonable to 
Beeume that the transverse equilibrium equation of Walters and Jones 
must be modified by inclusion of the term 

M9 My 


oo | ree (5.9) 


which is generated by the ~~ term. 
R2 
Hence, the modified transverse percuvogulenecquatton, becomes 


N | Wee -+- Z vee + 
Ae . + cotdw 4 csc* ow a9) M 


o> >1O9 





+ "2eEGEegM: eee Meee csc7 om! 


d,9 6 6,90 


+ 2cotécscdomM' + 2cscomM! cotpM! , = 0 (5.10a) 


66,0 a) 


6,60 — 


The in-plane perturbation equations, (5.6 b-c) are rewritten 


in terms of w as 


2 2 a S) Zar? ee 
5 csc“~odcotow ay + 5 csc“ ow $00 > 0 (5. 10m) 
and 


e 5 akan 
cotow', + (Go 4cot*$)w,, 


one) 


“+ C867 OW" aa + 5 escécotéw, 544 = 0 (5. 1Ge) 

To remain consistent, the preceding alteration is only valid 
if equations (5.10 b-c) are either equal to zero or nearly zero. The 
equations can be nearly zero since the transverse equation (5.10a) is 
considered the more dominant equation. 
Meee Linear Visco-elastic Case 

The dominant transverse displacement rate is achieved by sub- 
stituting the linear constitutive equation (3.3a) into the equilibrium 


equation (5.6a) 





We (5 cb) 


Utilizing the following Legendre identities 


lero + cotorn,, + (csc*o + An) Pn 0 (5. 12a) 


II 


o> > 


and 


Pn, + cot$Pn, + (m@csc* + An) Pn @) C52) 


It 


o> > 


where 


hr ee) 


=35= 





The perturbation equation (5.10a) can be recast to 


yee 2 (v2 en) 0 (G1) 
where 


v2() = (),,, + cote( ),, + esc*( ) cae (5.13b) 


> $ 


and 


e 


3aqaw 


2€ 





Substituting the assumed solution for the perturbation trans- 
verse displacement (5.1) uncouples the amplification function in the 


form of ordinary differential equations of constant coefficients 


e e 
?AnBy + oe An) (2 - An) Bo, @) (52144) 


0B tale’ — in) (2 eee 0 = (52a) 
mn 2 mn 
Equation (5.14b) is identical to equation (5.14a) if mis understood to 
include the case m equals zero. 
Sis POV Leng Po 1S “cons tan & 


Ant 


B = D exp 3 
(+ —- An) (2 - An) 


ae Bes (Seals) 


where a is an arbitrary constant dependent on initial imperfections. 
IS : ; ; : 
When x 1s defined as %t a non-dimensional plot, Figure 9, can 
be constructed to exhibit that the critical mode is n equal to two. 


Lt sene Lunction 
2 


a pi, (cos@) cosind (5 2G) 


2m ie) 


1S calculated for a number of points of various $6 and 6 values using 


computer program A of Appendix (A), the two dimensional contour surface 


Bee tg 





of Figure 10 is generated. On the contour map, positive values in- 
dicate an outward elevation of the mid-surface from the mean radius, and 
a negative value indicates an inward depression of the mid-surface from 
the mean radius, due to the perturbation solution. Since the elliptical 
mode has been proven to be the most critical mode, it will be henceforth 
assumed that the initial imperfections are of the elliptical shape. 
meee Non-linear Visco-elastic Case 

Substituting the constitutive equation (3-8a) into (5.6a) 


produces a dominant transverse displacement rate of 


. a25 ox (oD) 


v= (Se) 


In a method similar to the linear case (eqn. 5.12a), the 


perturbation transverse equilibrium equation is transformed to 


AV2w' — (V2 4+ =) (v2 Does (5.18a) 


with : 
astw (5.18b) 


— 
Hl 

i 
No] uw 


When the curvature change equations (5.4 d-f) are combined with consti- 

tutive equations (3.8 d-f) and then substituted into equation (5.10). 
Using the proper Legendre identities (5.13 a-b), and assumed 

Solution for w' (5.1), the wncoupled expressions for the amplification 


function become 


iI 


Vie ems On = nuie 0 (5 ao) 
on 2 on 


and 


0 (5.19b) 


ci reieeean mene Cfo be args 
mn 2 mn 
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Maus. previdirg De is constant 


ANT. 
Bo =D. xO eee ae (5.20) 
(— = kn) (2 = An) 
2 
where 
De is an arbitrary constant dependent on 
initial imperfections, 
and 


La nate ; 

The critical mode, as in the linear example, is n equal to two. 

As with the infinitely long cylinder, the assumption that the 
effective stress Hoeene constant must be modified. mPHactiae stress is 
assumed to be piecewise constant. Time is divided into a finite number 
of segments within which the effective stress is assumed to remain con- 
stant. When the effective stress increases by a given percentage, the 
Segment is terminated, and a new segment with an increased value of 


mmenreccive stress is initiated. 


From equation (3.5b) the effective stress equals 


o* = — (o, - oie +o 4 +o.% + 61 C520) 


Jr d S) 6 


where, from equations (1.2 a-f and 3.8 a~-f), the stresses are expressed as 


ws = a + oe (5.21b) 
or uaealtars + oh (5221) 
“40 = Eoay: (5.2 1d) 
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From equations (5.7 a-f) and 5.8 a-c) expressions (5.21b - 5.21d) can 
be expressed as 


. = 6. (5.22a) 


where H is the shell thickness of the sphere, 





ot = © kyo x OOM ee, A = em) (5.22b) 
rewritten in terms of w 
w' 
Dee leans ° ‘o 
t ee ele Ses ' 
O A 3 2 fo} CG w' + Ww $6 35 5 coté 
w', 





2 csc*$) (5.22c) 


which can be further simplified using Legendre identities (5.12 a-b) to 


ee ae (5.224) 
d 3R2 o> 


in a Similar manner 


_KH |, (1-™) 3 


' ae sie - . ary er Z 
O 9 ; (5 w' + Wa coe WwW rggcse db 
R 
a wit (Sneeze) 
2 o> 
which is Simplified to 
KH = ° e 
Oo" = ae M) Cog + Mra a (50222) 
3R2 
Finally the shear stress is written as 
K (1-M) - 
i = — o* (5222) 
¢o 3 Yh2 : 


Peoeenm, by substituting equation (5.8c), is 


2x Ho * 44M) 
3 4 = rect oon -~w', 


csc) (5.225) 
8 
3R2 ? 
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Since the perturbated strains depend only on curvature change 
rates (5.8 a-e), which are all expressed as functions of w', it is 
postulated that the spatial point on the sphere which exhibits the maxi- 
mum perturbated transverse deflection rate will be the point of maximum 
strains, thus stresses. To locate the point of maximum curvature rate (w") 
one can utilize a precedure to study maxima and minima of functions of 


several variables outlined by Hilderbrand (1962). 


(a) Evaluate Dp and oe and then set equal to Zero. 
(b) Evaluate BT, : ae and ae 
Ciel f DO I ae es Wag , and Wa > O ,.a relative 
minimum exists. 
(ii) if BURN ae > W' 49 , and a <0 , a relative 


mMoaxlMnum exises. 


From equation (5.1) the critical value of w'(n = 2) is 


wi =B (t) (P (cosé) + cos@P! (cosé) 
2m 2 zZ 
4 cos20P* (cos})) (5.23a) 
with 
PB (cos®) = > (3cos"6 + 1) (5.23b) 
GP» (cos®) 
1 _ : andes oS : 
P 5 Soeee sing eee 3singcosd (32 Se) 
a2P5 (cosd) 
pe (cos¢) = sin*¢ = 3sin2$ : (S7223¢) 
é a(cos$) 2 
anus, 
w', = -—3cososing + 3cos6(cos7¢ - sin“¢) 
+ 6cos20sindcosé (5.23e) 
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Wing = -3sin8sindcosd - 6sin20sin7¢ (52236) 


Both of these expressions (5.23 e-f) are zero when ¢ and 0 equal 


m/z 
Evaluating the second derivatives 
Wry, = 3(sin2¢ - cos7¢) - 12cos@singcos¢ 
+ 6cos26 (cos7$ - sin7¢) (5-234) 
wo = -3cos8sing¢cos¢ - 12cos20sin74 (5.23) 
and 
ee = +3sin0(cos7¢ - sin29) _ 12sin20cosdsing (52.230) 


at ¢ = 71/2 and 8 = 7/2 one finds 


w' a (5.244) 
“oo 
Ww a ae (5.24b) 
W AY: = =3 (52 24e) 
oa: cae - ee . e e * 
Thus, w aay ay: W 46 and a relative minimum exists at 6 


Poual to 7/2 and § equal to 7/2. Noting, however, that the deflections 
rates herein are of a negative sense (B is always negative), the point 
¢= 7/2, 9 = 71/2 is actually a relative maximum. Further, this conclusion 
isin agreement with the maximum of the deflection field illustrated in 
Figure 10, which was calculated numerically. 

Substituting values (5.24 a-c) into equations (5.22 d, £f, h) 
yields 


i aan e 
5) a J 2B KH GMs 


4 2 le (5: 25a) 
R 
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1. JRE ee (5.25b) 
5 62 (B2 2m 
<* eeee, SH oa mene (5.25) 
oe) RZ 2m 
From equation (5.15) é 
y) be te 
B' (t) = — — a 
a ) ere oo 
p ip 
S Reece ia te) S 
=> = ira p(s ae 
5 on 7 O exp (>) (5.254) 
the perturbated stresses simplify to 
Se a ‘s 
Gl ee es ae See 5 
4 a (pO R* a> exe ( e) 
p i : 
Oe = eee S 
Ce ee aes —) = -7o' 2 
6 7 = R x CP (> 0 ‘4 (5 = 2 Sit) 
Te = 3 (2) r* ae ) =.<46: (5225) 
$0 H k PG $ eee 
where 
ee maximum deviation from radius of imperfections 


radius 


Hence equation (5.21la) becomes 


ge =e |(80" 7 + (64 + ot 4)? + (0, - 7o',)4 


V2 : : 
213s 
+ Oe e ry? | (5.204) 
Simplified to 
at dio = 
* = —1|20,% - 120.0' 2000! 2]? : 5 .26b 
O al ae ou $ + O 5 | ( ob) 


Substituting expressions (5.22a and 5.25e), the effective 


Seress equals 


; 
1 PR E 9 OR 
ok = S/S = (SR exp( 
H |2 2 ‘H 3 
/2 “ 
: c 
+ 112.5 (2) ?R*exp? sy] egos 


The initial effective stress = = O) equals 
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2) i. 
of EPB(L- 2B ws 112.58) 2Rx2|"  (8.26a) 
0 /> lis) De eas H 
Defining 
* 
v= sae (5772 Ja) 
O 
0 
F.-= het (522 7) 
S 0 
and 
A = initial value of A (eqn. 5.18b) 
0 
the effective stress can be recast as 
M-1 
Ee aie ou ae tea pot 
o* = ———]|->— - = (—) R*v exp 
Jo a Z 2 H 5 
; z yin 5 
+ 112 (2) 2px2y 0?) exp” = “as (52270) 


AS an example, 


Wi those Re ieee Cecual to. 10, .005, 3, 
0 


1.1 respectively, the growth of the amplification function is as shown 


miragure J). Note that 


until the last ten percent of the lifetime of the sphere. 


the amplification function stays relatively small 


As expected, 


the time segment duration decreases as time increases. 


It should be further noted that the in-plane equilibrium equations 


(5.10 b & c) are equal to zero at the point of maximum stress and trans- 


verse deflection (9 = 71/2, 9 = 


se) ae 


AA 





Ge ~DESCUSS LON 


Pee infinitely Long Cylinder 
mei Comparison to Past Results 

The most convenient basis of comparison of the analysis herein 
with past investigations is the predicted critical time of creep buckling. 
The time behavior of the initial imperfection growth is such that it will 
become large and grow rapidly only for time near to the critical time. As 
an aid for the reader, Appendix B has been provided to transform the nota- 
tion of past analyses into the notation used herein. 

Hoff et al (1959) predict a critical time of 

1oa( + 3 /(2a%2)?) 


ae ea ee ee (627) 
bt a, (Mr2) 


p 
2han ‘h 

It is at this critical time at which the amplification of the 
initial imperfections becomes indefinitely large. The herein results have 
been compared to Hoff in three ways by varying the values of a/h, a*, and 
Min Figures 12, 13, 14 respectively. The comparison reveals favorable 
agreement for smaller a/h ratios (a/h < 25) with a* and M equal to .005 
and 3 respectively; and small imperfections (a* < .003) with a/h equal to 
lO and M equal to 3. For large a/h ratios, the Hoff estimate of buckling 
time is approximately an order of magnitude greater than the presented 
analysis. The Hoff solution is less sensitive to larger, gross initial 
imperfections, predicting a collapse time of five times greater than the 
een analysis when gross imperfections are assumed. 


The plot of varying creep exponent, M, (Figure 14) shows excellent 


agreement for M equal to one and seven and an estimate five times less than 
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Mert for M equal to three and five. In Figure 14, a* and a/h equal 
.005 and 10 respectively. 

The actual time history of the two amplification curves are 
hearly identical, as exhibited in Figure 15. The amplification of im- 
perfections remains very small until the last ten percent of the life- 
time of the cylinder. 

Tne discrepancy in predicted critical time of the two analyses 
do not seem too serious in the light of the fact that creep buckling 
times are of the order of 10° hours, so that excellent agreement is ob- 
served for the first 10+ hours. 

Berman et al (1974) presented the particular case of stainless 
Steel 304, which has the material properties at 1050°F listed in Ap- 
pendix C. One finds a great overestimation of Berman when compared to 
both the analyses of Hoff and herein. Berman predicts a buckling time 
Berecylinder of a/h equal to 10 and a* equal to .05 of 1.743 x Lo“ p88 
mate Hoff predicts 3.6852 x 10° po ae and the herein investigation 
estimates 1.826 x ioe pane. Both the Hoff and herein analysis 
predict a buckling time of approximately four orders of magnitude sooner 
than Berman. The author feels that the calculated time of Berman is in 
possible error, Since he uses an out-of-roundness ratio of .05, a very 
gxoSsS out-of-roundness for any realistic engineering design. Further 
Berman fails to update the effective stress with time. The herein 
analysis does predict a buckling time of 8.3 x 10 reap ee ripe fel aS 
005, (an order of magnitude less) and the effective stress is not up- 


dated with increasing time. 
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The numerical analysis of Samuelson (1970) deals with the creep 
response of finite length cylinders subject to both non-uniform and 
uniform loads. Thus, higher harmonics come into an influencial role in 
the creep process. Further the transverse deflection is dependent on 


fe Jongitudinal, x, coordinate 
w' = be A (t) cosn@sin ues (622) 
=O n L 


Samuelson predicts a buckling time of approximately 100 
seconds for a cylinder with a/h, L/a equal to 27 and 2 respectively. 
A similar buckling time is found for the same cylinder by the herein 
analysis if M equals 3, a nominal a* .002, is assumed and < is 

a 

assumed to be of the order 10°. ‘Thus, it is possible to state that the 
herein analysis predicts a buckling time of the same order of magnitude 
as Samuelson. 
6.1.2 Elastic Response 

The elastic response of the cylinder can be obtained by sub- 
stituting the elastic part of the constitutive relation (3.2a) into the 


transverse equilibrium equation (4.6a) to yield 


2 ae 
a p, (1 U~) 


ee = = (6. 3a) 


The elastic transverse deflection from Roark (1965) is 


a°p, V 
WwW = (1 - —) (6. 3b) 


e kh 2 
Roark 





If Poisson's ratio, Vv, equals .3, the ratio of the herein 


elastic displacement to that of Roark is 1.07. 
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6.1.3 Parametric Study 

Due to the dependence of the time history of the transverse 
deflection on a multitude of independent variables, radius to thickness 
ratio (a/h), out-of-roundness ratio (a*), material constants (M, kK) and 
loading (po); an indefinitely large number of parametric curves are 
required to adequately explore the creep response of cylinders under 
extermal pressure. This task is greatly simplified by the computing 
of transverse displacement (4.18a) using computer program B of Ap- 
pendix A. 

Figures 16, 17, and 18 are plots of non-dimensional critical 

Bees ; 

time, ee lu versus a/h for creep exponent M equal to 3, 5, 7 re- 
Spectively. In each case as a/h ratio for a given a* decreases, the 
non-dimensional critical time drastically increases. Further, as a* 
increases for a given a/h ratio, the non-dimensional critical time de- 
creases. Of the two parameters, a/h and a*, a/h is the more dominant 
Seeraole. That is, critical time is more sensitive to a/h than a*. The 
Same conclusion was reached by Hoff (1959). 

Figure 19 illustrates the effect of varying the initiation 
point of a new time segment of constant effective stress. If for example, 


t gt 


Sc rary 
ee ea S 


ieeeey ni, a* equal 3, 10, .005 respectively, the value K 


3.5 x 10 © for time segments which are initiated when the effective 
Stress (4.19) increases by fifty percent () = 1.5) compared to a value 
of 2.2 x 10 © when the time segmentsare initiated when the effective 
Beress increases by five percent (jf = 1.05), a ratio of 1.6 difference. 
Froi1 a practical ocean engineering point of view, a comparison 


of the two most common creep prone materials, aluminum and titanium, used 
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in the ocean is in order. The alloys used in the comparison are Aluminum 
6061-T6 and Titanium 6Al-2cb-1Ta-.08Mo. The room temperature creep 
constants of these materials are listed in Appendix C. The titanium con- 
Stants where experimentally found by Chu (1972), while the aluminum con- 
stants by the author. 

The first comparison is non-dimensional critical time verses 
the initial elastic hoop stress, expressed as a percentage of yield stress. 
The yield stresses of Titanium 621/.08 and Aluminum 6061-T6 are 120 and 
40 ksi respectively. The two materials are compared with two a/h ratios, 
10 and 50; and two a* ratios, .Ol1, and .OO1. A number of conclusions 
can be drawn. The lifetime increases for both materials, env oxiinacee 
mye times for a/h equal to 10 and forty-five times for a/h equal to 50, 
with an order of magnitude decrease in initial out-of-roundness ratio 
(a*). Hence, it appears to be more important that cylinders of large 
a/h ratios to be constructed to closer tolerances than cylinders of a 
smaller a/hn ratio. Titanium is the better performer at low stress levels, 
but is surpassed by aluminum at higher stress levels. Imperfections play 
a more important role at higher stress levels. Finally, for a stress 
meyel less than forty percent of yield, there is a dramatic increase in 
lifetime. 

mie Ioeivuence the varying a/n at a Gavem seress level Sixty 
percent of yield, has on the non-dimensional critical time is portrayed 
in Figure 21. Again, a crossover occurs. For small a/h ratios (a/h < 20) 
titanium is the optimum material, while for larger a/h ratios aluminum 


exceeds titanium. 
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Taking another tack, if a* is varyed for a constant a/h (10) 
and a constant stress level, fifty percent of yield, Figure 22 is gen- 
erated, a crossover of performance is observed with titanium being 
Surpassed by aluminum for grosser tolerances (a* > .007). 

The most pertinent collation is the comparison of the two 
materials at equal pressure levels. As an example, with Pi! ax, a/hi 
equal to 2000 psi, .001 and 10 respectively titanium out performs 
aluminum by four orders of magnitude, (Figure 23). Titanium, with its 


Superior strength, is able to operate at a much lower percentage of 


yield stress. 
pee Complete Sphere 
fee-L Elastic Response 

By substituting the elastic part of the constitutive equation 
(3.2a) into the transverse equilibrium equation (5.6a), an expression 
for the elastic response is developed 


a*p (1 - 0%) | 
v=o (6. 4a) 


@ Ben 
The elastic deflection given by Roark (1965) is 
a*p (1 - v) 
ue = SER (6.4b) 
Roark 

A ratio of the elastic deflection to that of Roarkis .86 of v is 
assumed to be .3. 
6.2.2 Parametric Study 

As with the infinitely long cylinder, the dependence if the 
time history of the transverse deflection on a large number of indepen- 


dent variables requires the generation of a large assemblage of para- 


metric curves to adequately investigate creep response of the sphere. 





: 





Again, the task is simplified by numerically computing equation (5.1) 
using computer program C of pee eet. 
bs Figures 24, 25, and 26 are graphs of non-dimensional time 
= versus R/H for various creep exponents, (M equal to 3, 5, 7). 
The non-dimensional critical time decreases rapidly with both increas- 
ing R/H and R*, being most sensitive to R/H. 
The non-dimensional critical time of the sphere is fairly 
non-sensitive to the duration of the time segments, as illustrated in 


Figure 27. If the updating of effective stress (5.26a) is carried out 


when effective stress is increased by fifty percent (i = 1.5) the 


ceeds - 
malue of — a «As lees 10s. If, however, the Etec is updated 
: iE 
when it has increased by five percent, ( = 1.05) aE cise ales eS 


x 10 °, an increase of flit tyepenrcent- 

When the non-dimensional critical time is calculated for 
various initial elastic loadings, expressed as percentage of yield 
stress, for the two alloys, Aluminum 6061-T6 and Titanium 621/.08, 
Figure 28 is manifested. Initial imperfections are quite influencial 
on the critical time, causing a two order of magnitude increase for 
a/h equal to 10 and six orders of magnitude increase for R/H equal 
to 50 with an order of magnitude decrease in initial imperfections. 
Thus, revealing that for spheres of large R/H ratios, close tolerances 
are a prerequisite. Except for the case of high R/H ratio and gross im- 
perfections, titanium is the superior material, producing a longer 
lifetime. Rapid increase in lifetime is achieved if the stress level 


of the sphere is less than sixty percent of the yield stress. 


=o 





When R/H is varied for a given stress ie oT percent of 
yield, titanium, regardless of the R/H ratio, out-performs aluminum. 
This superiority is even = obvious, aS demonstrated in Figure 29, 
for substantial R/H ratios (R/H > 30). 

If for a constant R/H ratio and a loading of fifty percent 
of yield, the magnitude of initial imperfections is allowed to veer, 
Figure 30 is developed. Again, titanium is always the superior per- 
former, especially with small initial imperfections. That is, one 
has much more to gain by machining a titanium sphere to exacting 
tolerances than an aluminum sphere. 

At the same pressure level, the titanium is vastly superior 
than the aluminum, having a lifetime of approximately 7 orders of 
Magnitude longer with P., equaling 2000 psi and R* equaling .001. The 
obvious reason being that a titanium sphere will operate at one-fourth 
the percentage of its yield stress at a given ambient pressure. 

6.3 Comparison of the Creep Performance of the Sphere and Cylinder 

If the critical time of a cylinder and a sphere of the same 
material at equal pressure levels are compared by plotting the ratio of 
eeeeeal time of the sphere to the critical time of the cylinder for 
various a/h ratios, initial imperfections, and creep exponents, (Figures 31, 
32, and 33) it is discovered that regardless of the initial conditions, 
the cylinder has the inferior creep resistant properties. This point is 
especially manifested in the case of titanium, M equals 7. The oscil- 
ieerons observed in Figures 3), 32, and 33 result from the variation of 


the ratio of the effective stresses with various a/h and a* combinations. 





The sphere 1S more sensitive to variations in initial im- 
perfections, which is to be Sees Since the cylinder is a develop- 
able surface and is less Sees than the sphere. 

The sphere has excellent creep characteristics for stress 
levels less than sixty percent of yield stress. The cylinder, on the 
other hand, the stress level must be less than forty percent of yield 
Stress for favorable creep properties. 

Finally titanium, because of its high yield stress is in 
every instance the best material to be used in ocean engineering ap- 
mPlucGeti1ons of both the sphere and the cylinder. 

6.4 Design Criterion 

The extensive analysis present suggests an assemblage of 
Simple design guidelines to assist the engineer in developing a cost- 
effective, creep resistant design. In most ocean engineering shell 
designs, the major constraint in design is considered to be the minimi- 
@Zation of the weight to displacement ratio of the structure, while 
maintaining a maximum stress level of approximately fifty to sixty per- 
cent of yield. It is this constraint which dictates the use of light- 
weight, but creep-prone materials iyechew rl rsesp race. (remeunazely 
this criterion is satisfied if titanium is used rather than aluminum in 
pressure vessels, even though titanium is 1.5 times as dense as aluminum. 
The reason being, the aluminum sheil would have to be four times as 
thick to achieve a stress field of equal percentages of yield stress. 
Hence, titanium will always give a longer lifetime. 

While the sphere furnishes a vastly Superior creep response 


mmo the cylinder, it is more expensive to construct and machine. Thus, 





the engineer must weigh the merits of the costs of materials, con- 
GSeem@etion and machining of both an appropriate sphere or cylinder for 
his particular situation (pressure level and time constraints) before 
he develops a final design. 

Regardless, the shape of the amplification curves for both 
the shell shapes suggest a creep design limit. For approximately the 
first seventy percent of the lifetime imperfections remain small and 
grow in a linear manner. Thus, an arbitrary design limit can be estab- 
lished which states that as long as the imperfection growth becomes 
linear, the design will not creep buckle. The in-service monitoring 
of the creep deformation process is easily accomplished with a three 
element strain gage rossette placed at the point of greatest imperfec- 
tion. When the effective stress approaches the magnitude of the yield 
stress, the active use of the structure should be terminated. 

6.5 Example Calculations 
6.5.1 DSRV ALVIN 

The research submersible, ALVIN, owned by the Office of 
Naval Research and operated by Woods Hole Oceanographic Institution, 
has a spherical titanium 621/.08 alloy pressure hull. The 40.005 inch 
ragius hull has a a/h ratio of 20 and a maximum out-of-roundness ratio 
(a*) of 0.002. The maximum operating depth is 12,000 feet (5280 psi) 
and a collapse depth of 18,000 feet (7920 psi). 

The secondary creep law for Ti 6Al1-2cb-1Ta-0.8Mo alloy at 
room temperature has been empirically measured by Chu (1972) (Appendix C). 


While the measurements were done in tension, it is believed that the 





compressive creep behavior is similar to that of tension. The law is 
expressed as 


o/ Ke (Gee) 


where 


7 
ee 2 eyes so) CE ie (6.2b) 


in 

With values of a/h, a*, M equal to 20, .002, and 7 respectively, 
the amplification curve of Figure 34 is manifested. If subjected to 
the maximum operating pressure, buckling occurs at a time equal to 
9870 hours. If the submersible is operated at only half its operational 
depth it would not buckle until 1.85 x 10© hours had past. 
6.5.2 SOFAR Floats 

The Woods Hole Oceanographic Institution has employed a long 
cylinder which when released at sea sinks to a depth at which it becomes 
neutrally buoyant. The current measuring floats, first designed by 
Swallow (1963), will follow a water mass of a constant density due to 
the neutrally compressible, buoyant properties. Placed at a depth of 
approximately 5000 feet, (the axis of the SOFAR Channel, a horizontal 
channel, in which the convex shape of the sound velocity-depth profile 
allows sound to be transmitted with minimal attenuation over thousands 
of miles), the floats emit an acoustic tone at a given frequency. 
The float is then spatially fixed using acoustic triangulation from 
land based listening stations. The floats have been successful in 
tracking eddies in the WIOI MODE program. Unfortunately, after two to 
three months, these floats have been found to sink an additional 230 to 


960 feet, with accelerating sinking velocity, (an apparent creep buckling 





phenomenon). Mavor (1974) computed a linear strain rate of 0.9 win/in/hr 
from the change of displacement and an assumed constant in situ seawater 


density gradient. 


Ap “6 
=== = 6. 
DAL 1 o7 x 10: ert (G23) 
where 
0 = density 
Z = depth 


The floats are standard ALCOA extruded 6061-T6 aluminum 
cylindrical stock with a mean radius of 5.625 inches, a a/h ratio of 
7.5, an approximate out-of-roundness ratio of 0.005 and a length of 
17 feet. They operate at a mean depth of 5000 feet (2200 psi). The 
6061-T6 aluminum follows a cubic creep law aS experimentally found by 


the author (Appendix C) 


o3 = KE (6.4a) 
where 


4 x 10!9 (psi) 2hr. (6. 4b) 


7 
il 


Figure 35 displays the predicted creep response with a/h, 
fee equal to 7.5, 0.005, and 3, respectively. 

Since the problem is concerned with only the first three 
months of the lifetime, the perturbation transverse displacement can be 


considered linear. From Figure 35 


a SASS unlis, Jee (556 


From equation (4.1), the perturbation hoop strain rate is 


eae 2hw' 
) as: 2 
ac 





= 5.3947 x 10 8 in/in/nr. GS) 
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From equation (4.7b), the dominant hoop strain rate is 


e,.=~ (6.5c) 
ee a 
Computing w from equation (4.8) yields a dominant hoop 
strain rate of 
By = 9.45 x 10 ” in/in/hr. (6.5d) 
Hence, the total predicted linear strain rate 1s 
@p = ey + e', = -998 win/in/hr, (6.5e) 


ten percent greater than the estimate of Mavor. 
This calculation shows that, in fact, creep deformation is 


most probably the cause of the unexpected sinking of the floats. 





7. SUMMARY AND CONCLUSIONS 


The theoretical investigation which has been undertaken into 
the creep instability of the infinitely long cylindrical and complete 
spherical thin shells has satisfied the original thesis objective of 
developing a relatively simple analytical method to predict the creep 
deformation process of these shells. It was found that in both shells 
the critical buckling mode is the second elliptical mode. The roles of 
shell radius to thickness ratio and initial out-of-roundness have been 
explored. The relative superiority of the spherical shell over the 
cylindrical shell has been discussed. 

The investigation demonstrated that the single most important 
factor in creep design is designing the structure in a manner which 
Maintains the stress field at the lowest possible level, thus extending 
the structure's lifetime. This one factor endorses the use of high 
strength titanium and the spherical shell in creep design. 

In conclusion, while creep buckling can be a problem in long 
mem ocean Structures, the nature of the time history of imperfection 
growth, almost negligible until the last ten percent of the lifetime, 
allows the ocean engineer to develop a safe design resistant to creep 


Maexling with little difficulty. 





8. RECOMMENDATIONS 


Two essential advances in creep design are required in the 
future. First, the establishment of accurate creep constants at 
ambient temperatures for common ocean engineering materials. Second, 
a more advanced theoretical investigation which encompasses the 
effects of boundary conditions of such shells as finite length Saline - 


and partial spherical shells (end caps). 
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APPENDIX VA) = COMPUTER UTELI ZATION 


feeecomputer Program A - Perturbation contour map of the complete 
sphere. 


The following FORTRAN program was developed to compute 
wi /Bo for a discrete number of points, (Ss by 5° grid) in a given 
quadrant. One needs only to input the initial values of $ and 90 


of a given quadrant to receive the x-y position and value of w'/B, 


of the before mentioned grid. 


x = sindcosé (1) 
y = cosé (2) 
w' zZ m 
— se P 5 (cos$) cosmé (3) 
2m m=O 


Peet URBATION CONTOUR MAP OF THE COMPLETE SPHERE 
DIMENSION T(19), P(19), W(19), E(19), F(19), P2(3) 
DO 51 M=l1, 4 
READ(8,200) U,V 
200 FORMAT (2F10.6) 
bO 52 N=1, 19 
A=N-1 
T (N)=.0174e5.0¢eA+U 
DO 53 I=1, 19 
B=I-1 
P (I) =.0174e5.0eB+V 
S=SIN(P«1) } 


C=COS (P(T) ) 


= IOS 





P2 (1) =.5¢(3.00Ceo2.0-1.0) 
P2 (2) =3.09S¢c 
P2 (3) =3.00See2.0 
wtr) =020 
DO. 4a — lms 
R=J-1 
D=COS (T (N)e R) 
WZ=DeP2 (J) 
54 W(I)=W(I) +wWz 
E (I) =SeCos (T(N) ) 
F(I)=C 
iremme (5,202) E(1), F(t), W(t) 
Zee rORMAT('', 3F15.5) 
Eee CONTINUE 
ee CONTINUE 
51 CONTINUE 


END 


Notation Of Program A: 


T 8 
P p 
BA {CIL) P, (cos) 
E2K2) P*, (cos$) 
Zz 
P23} P 5 (cos) 
W transverse displacement/amplification function 
E | x 


Oa 





PB y 
U initial 6 value 
V initial ¢ value 


ieee computer Program B - Growth of the amplification function of the 
intinicely wong cylinder. 


This program was devised so that one inputs a/h, a*, M and 
uy, and equation (4.19) is computed. The output consists of non- 
dimensional time oe normalized amplification function as oe and 


a second non-dimensional eiMe wibcnews.a betters 1naieator of reall 
tP 


k 





time 


INFINITELY LONG CYLINDER 
MS=3 
ie 15 
pool M=1, MS 
READ(8,100) B,D,E,XM 
100 FORMAT (4F 10.6) 
T=0.0 
Y=1.0 
WRITE(5,101) B,D,E,XM 
omnonvct@e |! Sx vw’ FiO GA =", FlOso, lOx a Phi"  FlOcG, 10x72 M 
l= ee hOeG. // ks eet Me 20. AME EN Ore Real ot ME 4//) 
bere d=1,0S 
Pe S(inGl. ))- GOrrear 
=I 10 
GO TO 8 
(ai ex 


F=Yeoe(XM-1.0) 


= U0) Bi 





8 S=EXP (TeF/3.0) 
A=E @(1.0+6.00eEeBeDeEXP (TeF/3.0) ) 
T=3.0/Fe°ALOG ((A-1.0) /(BeDeFe6.0) ) 
T1=T/(9.0 Bee(XM+2.0)¢(.8660(1.+6. B D) )e e(XM-1.)) 
WREEE (S03) T,A,;T1 

HOG. PORMAT (185). 7xX 610.3 ,1/% ,F1l0. 3, 15%, 82075) 
27 
1 CONTINUE 


END 


The number of input statements is controlled by the variable 
feeewiie the number of cycles for a given set of inputs is controlled 


ly the variable JS. 


Notation of Program B: 





T T 
Cc 
A A /C 
2 2 
ue 
ard: i 
oe p 
XM M 
B a/h 
D a* 
E 
is 


(iii) Computer Program C - Growth of the normalized spherical amplifi- 
Cation LUuNeCtIOn wilhetime. 
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This program is structured exactly as the cylindrical use. 
Equation (5.1 ) is computed from the inputs R/H, R* M and ae 
The output received is non-dimensional time ode amplification, fumetion 
(BL /D5_) and a second non-dimensional time which is a better indicator 


(eae 
of real time = - Again, the number of input statements and number 





of cycles for a given set of input is controlled by the variables MS, 


JS respectively. 


COMPLETE. SPHERE 
MS=8 
JS=50 
DO 1 M=1,MS 
READ (8,100) B,D,E,XM 
100 FORMAT (4F10.6) 
m=O.0 
eae) 
WRITE (5,101) B,D,E,XM 
ore SOrmaAr(']' 5X, A/H=",F10.6,'A=',F10.6,10X,'PHI=" ,F10-6,10x,'M 
fete. 6 //, 5X7 TIME 720%, (AMP EN', OMe REAL TEMES) 
Domes —1 5 
PomeeGles|) GOVlo 7 
=. 
GO TO 8 
7 Y=ESx 
F=Y¥ @ o(XM-1.0) 


8 S=EXP (TeF/3.0) 
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R1=4.5eBeDeFres 
R2=112.5¢(BeDeFeS)e °2.0 
GaE 002 .0°(0.54+R1+R2)-0.5 
A=((-RL/S)+((R1/S)e 02.0+4.0°GeR2/Se 02.0)e 60.5) /(2.0¢R2/Se 2.0) 
T=3.0/eeaL0c (A) 
Tl=Te2./(9.eBee(XM+2.)0(1./2.¢ ©,50(.5-4. 5¢Bopt+112.5¢(BeD)e e2.)) 
0.5¢e(XM-1.) 
WRITE (5,103) T,A,T1 
ies FORMAM Ges P10,3,17X,F10.3,15xX,E20:.. 5) 
2 X=Y 
1 CONTINUE 


END 


Moeatiton Of Program C: 





T T 
S 
y p 
B 
E om/ 2m 
P 
Tl 7 
XM M 
B R/H 
D R* 
E p/w 
/ | 
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APPENDIX B - EQUIVALENT NOTATION FROM PAST INVESTIGATIONS 


femeott, N.J., Jahsman, W.E. and Nachbar, W. “A Study of Creep 


Collapse of a Long Circular Cylindrical Shell Under Uniform 


External Pressure", J. of Aero Space Science, Vol. 20, No. 10, 


1959. 
Notation of Paper Equivalent Notation 
OJ a* 
aR a 
= Qak(S 
x 3 a (7) 
meg 
== nk, ,PR, ,R,-2 eee, 2a. 0 
eas jee Oe eet pats Pas ill psa 
g ag) Mas ve a ae 
9 oh 2 
ikeie(UiGe = a =o 
1 9 8 a ao 
bin cae er log (1 +227 x 2) ee 
CE 2 2° 0 re 
A igual 
= K 
r 
n M 


h ; ; 
where 9 = 5 radius Ol soy baclon 


oo | 
H | 


average radius 


a = Shape factor 


~ 
It 


weighted shape factor 
meneeeerman, L., Chern, J.M. and Gupta, G.D. “A Parametric Study 
on Elastic-Plastic Creep Buckling of a Thin CyJandrical Shell’, 


B=) or Pressure Vessel Technology tugusrel ¢- 
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NGtateronmet Eaper Equivalent Notation 


ww 
ee, * 
oF, Da a 
a a 
h h 
gq , P 


z 


(iii) Samuelson, L.A. “Creep Buckling of a Cylindrical Shell Under 


Non-Uniform External Loads", Int. J. Solids Structures, Vol. 6, 


no70. 
Notation of Paper Equivalent Notation 
d 0 
B K 
R a 


The remainder of the notation iS consistent. 


sO Si 





APPENDIX Cae CREEP CONSTANTS OF COMMON 
OCEAN ENGINEERING MATERIALS 
(i) Titanium 621/.08 (Ti-6AL-1Ta-.8 Mo) 
(Oo = 120 ksi) 
y 
Chu (1972) found that the creep data of T1 6217-208 ateroom 


temperature can be represented as 


o  .13.890. 0.183 
€ = (03.55) = (1) 


where 


m 
ll 


creep strain, % 


stress, ksi 


0} 


1G 


time, hours 

L£ equation (1) is plotted at a given stress level, it is 
Giscovered that after an initial period of primary creep, the non- 
linear viscous behavior of secondary creep becomes apparent 

as example: for o = 103.55 ksi, € = .1071 pin/in/hr. 

By calculating the constant strain rates for a number of 
stress levels the following power law emerges 


o7 = 8.3 x 104% (2) 


where O 


stress, psi 


€ creep Strain rate, pin/in/hr 
(ii) Titanium 6-4 (6AL-4V) 
(Oma so, koa) 
Y 
Chu has found the emperical equation for titanium 6-4 at 


room temperature, but has not as yet published it. 


=o 





(iii) Aluminum 6061-T6 (wrought alloy - magnesium and silicon alloy 
solution treated and artificallyeacqed) 
(o st 40 ksi) 
y 
The author conducted a number of tests in which cylindrical 
columns of Aluminum 6061-T6 were loaded with a uniform load in compres- 
sion by a lever system. Deformation was measured using two uniaxial 
strain gages mounted on opposite sides of the column in a series con- 
figuration. The series configuration allowed bending effects to be 
internally cancelled out of the readings by the gages. 
The Aluminum 6061-T6 was found to abhor to a cubic law at 
room temperature 


o3 


4x 101% (3) 
where 


O = stress, psi 


creep Strain rate, pin/in/hr. 
All tests were conducted at high load levels, greater than 


65 percent of yield, at room temperature. The strain rates computed 


abe 
o (psi) € (yin/in/hr) 
34670 0.8504 
33750 0.8420 
31450 0.7768 
27615 0.0569 


me Stainless Steel 304 (Austentic) 


Be rman OES 7/25) expresses the secondary creep behavior of 


Stainless Steel 204 at 1050°F as 


=hLO= 





g2 6/89 2 apes e2 


where 


Q 
| 


stress, psi 


m 
H 


creep strain, pin/in/hr. 
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